Destruction of integer quantum Hall effect at strong disorder: A numerical study 
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A breakdown of integer quantum Hall effect (IQHE) at strong disorder is studied numerically in 
a lattice model. We find a generic sequence by which the integer quantum Hall plateaus disap- 
pear: higher IQHE plateaus always vanish earlier than lower ones. We show that extended levels 
between these plateaus do not float up in energy but keep merging together after the disappear- 
ance of plateaus, which eventually leads to a localization in the whole system. We also study this 
phenomenon in terms of topological properties, which provides a simple physical explanation. 
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It is an interesting question to ask how the integer 
quantum Hall effect vanishes in the strong disorder or 
weak magnetic field limit. Many years ago, Khmel'nitzkii 
and Laughlin |^ both argued that extended levels at 
the centers of Landau levels would not disappear dis- 
continuously nor merge together at strong disorder, but 
rather float up continuously to pass the Fermi level. This 
argument is a key scenario in the global phase diagram 
proposed by Kivelson, Lee and Zhang |3| for the quantum 
Hall effect. Recently a considerable number of exper- 
imental measurements [^-^ have been devoted to this 
issue, in which a limited floating-up feature [^[jsj is in- 
deed observed for lower extended levels. However, such a 
floating up seems not indefinite and experiment |^ sug- 
gests that higher extended levels eventually merge with 
the lowest one at weak magnetic fields. A direct transi- 
tion from a higher IQHE state to an insulator state has 
been also found. ||^,^ These are in contrary to a simple 
floating up picture and imply that an overall destruction 
of IQHE may be more delicate. Recent numerical study 
of finite-size localization length by Liu, Xie, and Niu ||] 
has also suggested no global floating up for extended lev- 
els, although no merge of extended levels was shown and 
actually how each of extended levels is destroyed by disor- 
der remains unclear. Thus the problem of disappearance 
of the IQHE needs to be reexamined in a more careful 
and global way. 

In this Letter, we present a numerical study of these 
issues based on a tight-binding model. A systematic de- 
struction of the integer quantum Hall plateaus is un- 
ambiguously revealed, which exhibits the following se- 
quence: higher plateaus are found to be destroyed first 
by strong disorder in a one-by-one order, and the lowest 
plateau is the last to be gone before the whole system 
becomes an insulator. In particular, the lowest extended 
level located below the last plateau is not floating up in 
energy, and even the critical exponent at such a delocal- 
ization point remains essentially the same as in the full 
IQHE case. A topological reason for such a destruction 
of the IQHE is given based on the Chern number which 
characterizes extended states | p^ : nonzero Chern num- 
bers with opposite sign are moving down from the band 



center to annihilate those in the lower-energy extended 
levels, and eventually a total annihilation of Chern num- 
bers leads to a global insulating phase at strong disorder. 
Furthermore, it is indeed found that a limited floating- 
up picture for the lower extended levels shows up when 
the electron filling number, instead of the Fermi energy 
itself, is flxed with the increase of disorders. At higher 
extended levels, we see a direct transition from the IQHE 
states to insulator states. 

We consider a tight-binding lattice model of noninter- 
acting electrons under a uniform magnetic fleld and dis- 
orders. The Hamiltonian is defined as follows: 
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Here ef is a fermionic creation operator, with < ij > 
referring to two nearest neighboring sites. A uniform 
magnetic fiux per plaquette is given a,s (f) = J2o ^ij — 
27r/M, where the summation runs over four links around 
a plaquette. We study the case in which the integer M 
is commensurate with the system width L. And Wi is 
a random potential with strength \wi\ < W/2, and the 
white noise limit is considered with no correlations among 
different sites for Wi. 

Fig. 1 shows an overall picture for the Hall con- 
ductance calculated by the Kubo formula, with a flux 
strength = 27r/8 at a 16 x 16 lattice size (only E < 
part is shown in the figure due to the antisymmetry of 
the Hall conductance on the two sides of the band cen- 
ter E = 0). In the weak disorder case {W = 1), there 
are three well-defined IQHE plateaus at cth = ve^/h 
[v — 1,2,3), corresponding to four Landau levels which 
are centered at the jumps of the Hall conductance at 
-E < 0. With disorder strength W varying from 1 to 6, 
we see a systematic destruction of these IQHE plateaus. 
At = 2, the third plateau (closest to the band cen- 
ter) begins to disappear. An increase of to 3 and 
4 will further result in the disappearance of the second 
IQHE plateau, while the lowest plateau still remains well 
defined. In the inset of Fig. 1, such the lowest IQHE 
plateau at = 4 is shown with sample sizes varying as 
8 X 8, 16 X 16 and 24 x 24. While the quantization of this 
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V = 1 IQHE plateau remains at different lattice sizes, 
the region of the transition towards zero is continuously 
narrowed with all the curves crossing at a fixed-point cor- 
responding to the Hall conductance ^e^/h, which should 
be extrapolated to a sharp step of jump in the thermo- 
dynamic limit. This resembles a typical scaling behavior 
found for the IQHE in the weak disorder case and the 
lattice-size-independent fixed-point represents the mobil- 
ity edge of the delocalization region. (Later it will be 
shown that even the critical exponent at this delocaliza- 
tion energy point remains essentially the same as in the 
full IQHE case.) On the other hand, the Hall conduc- 
tance in the higher energy region at W = 4 has no more 
IQHE structure and is seen to continuously decrease with 
the increase of the lattice size. A breakdown of the lowest 
IQHE plateau is eventually found at larger W^s {W — 6 
case is shown in Fig. 1). Here the Hall conductance in 
the whole energy region is much less than /h, which 
continuously decreases at larger sample sizes (from 8x8, 
16 X 16 to 24 X 24), indicating that it will scale to zero in 
the thermodynamic limit with all states being localized. 

To further inspect the critical behaviors shown in Fig. 
1, we use a different finite-size scaling method in terms 
of the Thouless number , [p^ This method has been pre- 
viously used by Ando [|13[ to study a similar problem, 
but much larger sample sizes (from 16 x 16 to 96 x 96) 
with more random configurations (200 — 30, 000) can be 
achieved here by using Lanczos diagonalization method, 
in order to accurately decide the critical points and ex- 
ponents. Fig. 2 shows localization lengths versus energy 
obtained by a finite-size scaling calculation. |Q Only the 
data for W = A and 5 are presented here where only the 
lowest IQHE plateau still remains. At = 4, Fig. 2 
shows that the localization length follows a scaling be- 
havior £^ '-^ \E — Eci\~^ with the exponent a; = 2.4 ± 0.1 
(which is the same as in the case when full IQHE plateaus 
are present at weak disorders jl^) on the two sides of 
the divergent point Ed — —3.40. Such a scaling law ver- 
ifies the existence of a delocalization fixed-point shown 
in the inset of Fig. 1 and also confirms that the lowest 
IQHE plateau at = 1 is robust as delocalization point 
Ed is well separated by a localization region from the 
higher-energy part. The localization length ^ exhibits 
an another dramatic increase which can be extrapolated 
as \E — Ec2\~*'^ when one approaches to Ec2 — —2.46 
from below. Such a second delocalization point Ec2 is 
more clearly shown at W — 5 case in Fig. 2. Above Ec2 
one finds a large but finite localization length which is 
rather flat up to the band center and whose magnitude 
is reduced at Vl^ = 5. A recent numerical calculation by 
using one-parameter finite-size scaling analysis also sug- 
gests such a region to be localized. Thus our results 
show that after losing the IQHE plateaus at higher en- 
ergies, original extended levels within this region merge 
and move down to Ec2 aXW = A and 5. Fig. 2 also shows 
that Ed is only slightly reduced (from —3.40 to —3.48) 



with essentially the same critical exponent x when W is 
increased from 4 to 5, whereas Ec2 is much more quickly 
moving down (from —2.46 to —2.96). Correspondingly 
the V — 1 IQHE plateau between Ed and i?c2 becomes 
narrower, and eventually these last two extended levels 
merge and disappear at a critical Wc- We find that the 
localization length becomes finite at = 6 in the whole 
regime, which is consistent with the conclusion drawn 
from the Hall conductance calculations, and places Wc at 
a range: 5 < Wc < 6 . At smaller disorder strengths, like 
= 3, we also see that more divergence points emerge 
in ^, in correspondence with more extended levels and 
additional IQHE plateaus shown in Fig. 1. 

A key thing to understand the above evolution of 
IQHE plateaus is the localization-delocalization transi- 
tion, and it is a well-known fact that an extended state 
can be characterized in terms of a topological quantity, 
namely, nonzero Chern number . And the boundary- 
condition-averaged Hall conductance is a summation of 
all the Chern numbers carried by states below Fermi sur- 
face: < aniEf) e'^/hY.e^<Ef C^"''- Here the Chern 
number C*'™^ is always an integer. Since its distribu- 
tion in the thermodynamic limit decides delocalization 
regions, one may define states with nonzero Chern inte- 
gers as extended states. Plots of the density of 
extended states p^xt are presented in Fig. 3a with the 
same flux strength as in Fig. 1. In the weak-disorder 
case [W — 1), well-defined peaks of pext are at centers of 
Landau-level bands, separated by localized regions rep- 
resented by plateaus in Fig. 1. Widths of these peaks 
will approach to zero in the thermodynamic limit, p^ ] 
Total Chern number for each of three lower-energy peaks 
is found to be exactly -1-1, which is the reason leading to 
three quantized Hall plateaus at +1, -1-2, and -1-3 in unit 
of /h shown in Fig. 1, when the Fermi surface is located 
between these peaks. The last peak closest to the band 
center in Fig. 3 carries a total Chern number —3, which 
guarantees that the Hall conductance in Fig. 1 falls back 
to zero beyond the third plateau when the Fermi energy 
approaches E — 0. This is a peculiar feature of a lat- 
tice model since when the whole band is half-filled, the 
Hall conductance has to be zero. With the decrease of 
the fiux strength and increase of the number of Landau 
levels or peaks of pcxt at the center of the Landau levels, 
one always sees that the last peak near the band center 
carries a large Chern number with opposite sign and a 
magnitude equal to the total sum of those at lower-energy 
peaks. 

In Fig. 3a, upper two peaks at W — 2 start to merge 
due to the disorder scattering. AtW = 3 they are pretty 
much mixed together. Because these two peaks carry 
total Chern numbers of -|-1 and —3, respectively, an an- 
nihilation of these opposite-sign Chern numbers occurs 
with the disappearance of the third plateau in Fig. 1, 
which leads to a substantial reduction of the magnitude 
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of the Hall conductance in that region. In fact, the over- 
all magnitude of the Hall conductance beyond the sec- 
ond plateau is less than the value of the second plateau 
at = 3, suggesting that the negative Chern numbers 
have already moved down from the original —3 peak and 
dominated this region. Notice that the lowest peak near 
—3.4 in Fig. 3a still remains separated from the rest 
spectrum at W = A. The inset in Fig. 3a shows such a 
peak becomes narrower and sharper with the increase of 
lattice sizes. This is consistent with the finite-size scal- 
ing result in Fig. 2 that a well-defined mobility edge still 
exists at Ed = —3.40. The inset of Fig. 3a also shows 
a second peak emerging near —2.5, which coincides with 
Ec2 in Fig. 2 and implies that negative Chern numbers 
would be located in this neighborhood at a large lattice 
size, representing a delocalization region beyond which 
the Hall conductance jumps back to zero from the v = 1 
plateau. Thus when the Fermi energy is well above Ec2, 
the system becomes an insulator with zero Hall conduc- 
tance. In other words, a direct transition from original 
1/ — 2 IQHE state to an insulator state is realized here, 
which is consistent with the experimental finding. 
Finally, at = 6, the negative Chern numbers reaches 
to Ed and annihilates the last -f 1 Chern number peak in 
Fig 3a, where we find a monotonic decrease of pe^t with 
larger lattice sizes which corresponds to a localization in 
the whole region as shown by previous analyses. 

We do not see any floating up of extended levels before 
they vanish. In fact, extended levels are all moving down 
at strong disorders except Ed which is essentially not 
changed until this last extended level is completely an- 
nihilated by negative Chern numbers moving down from 
the band center. Nevertheless, a limited floating-up pic- 
ture may still be seen in the case if the Landau-level 
filling number is fixed as shown in Fig. 3b. If we look at 
Pext within the lowest Landau level (with the occupation 
number ul < 1), the peak position does shift towards 
higher occupation number at a larger W . Since the peak 
of Pext will represent a mobility edge in the thermody- 
namic limit, this upward shifting implies that for a given 
filling number in the lowest Landau level, the extended 
states can float up to pass the Fermi surface when the 
disorder is strong enough, leading to an insulating tran- 
sition. Such a "floating up" is due to the reason that the 
number of localized states below the lowest extended level 
at Ed is increased at stronger disorder, so that when the 
electron filling number is fixed the Fermi level is actually 
moving down, or relatively, the extended level is shifting 
up. This effect is indeed consistent with experimental 
observations for lower extended levels. 

Hall conductances at weaker magnetic fields (thus with 
more Landau levels as well as IQHE plateaus) have been 
also calculated with flux strength cj) = 27r/M at M = 
11, 16 and 24. All of them exhibit the same generic 
features as shown in AI ~ 8 for the destruction of the 
IQHE at strong disorder. It suggests that at large M 



case (weak- field limit), lower extended levels are always 
fixed at lower energies until higher extended levels move 
down to merge with them after higher IQHE plateaus 
are destroyed in a one-by-one manner. In completing the 
present work, we also became aware of a recent Chern 
number calculation by Yang and Bhatt ||l^ at M — 3, 
by which they concluded that a "floating up" picture is 
correct. Such a floating up of the lowest-delocalization 
region seen by them is actually the same as shown in Fig. 
3b. But we point out that M = 3 is a very special case 
where only one plateau = 1) exists aA, E < 0, and thus 
one cannot find the sequence of the destruction of higher 
plateaus as well as the merge of higher delocalization 
levels, which are essential in our discussion of the global 
IQHE evolution at strong disorder. 

In conclusion, our numerical calculations unequivo- 
cally demonstrate that the integer quantum Hall effect in 
a lattice model is destroyed by strong disorder by the fol- 
lowing sequence; higher integer plateaus closer to band 
center are first to vanish when lower plateaus still re- 
main well defined. Once an IQHE plateau disappears, 
two delocalization levels separated by it will merge to- 
gether and move down. Since extended states at the band 
center carry a total Chern number with opposite sign as 
compared to those in the lower delocalization levels, the 
sequence of the disappearance of the IQHE plateaus ac- 
tually corresponds to that Chern numbers at the band 
center continuously move down towards the band edge 
and annihilate those Chern numbers with opposite sign 
in all the lower extended levels. No floating up in energy 
is found for extended levels during such a destruction 
procedure. However, a relative upwards shifting of lower 
delocalization levels is indeed seen if the electron filling 
number, not the Fermi energy, is fixed. Our results thus 
provide a consistent explanation of both "fioating up" as 
well as direct transitions from higher IQHE to insulator 
states in the experimental measurements. 
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Fig. 1. The Hall conductance (Jh as a function of 
energy E is plotted for different disorder strength W's, 
at a lattice size 16 x 16. The inset shows the evolution 
of the V ^ \ IQHE plateau aX W — A with lattice size 
varying form 8 x 8 (o), 16 x 16 (•), to 24 x 24 (+). 

Fig. 2. Localization length ^ obtained by a finite-size 
scaling based on the Thouless number is shown as a func- 
tion of energy E. Critical behaviors near divergent points 
are fitted (dotted curves) by scaling laws (see text). 

Fig. 3 (a) The density of extended states pext versus 
E. Disorder strength W^'s are the same as in Fig. 1 
with a lattice size 8x8 (about 400 random-potential 
configurations are used). The inset shows pext at = 
4 with lattice sizes varying as 8 x 8 (•), 16 x 16 (x), 
and 24 x 24 (*); (b) p^xt versus Landau-level occupation 
number til with the same parameters as in (a). 
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